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EXECUTIVE  SUMMARY 


'Conventional  linear  elastic  fracture  mechanics  cannot  explain  fracture  in 
concrete  and  mortar.  Unlike  traditiondi>4J2FM  materials  which  show  a  constant 
fracture  toughness,  concrete  and  mortar  tdtrghen  (show  more  resistance  to 
fracture)  as  cracks  propagate  until  an  instability  results  and  the  structure 
fails.  This  increase  in  toughness  is  thought  to  occur  because  of  many  possible 
complex  mechanisms  which  take  place  during  fracture  in  concrete  and  mortar. 
These  mechanisms  include:  formation  of  a  zone  of  microcracking  around  the  crack 
tip;  crack  bridging,  where  aggregate  particles  span  the  crack  and  act  as 
ligaments  to  hold  the  crack  together;  and  crack  deflection,  where  the  crack 
changes  direction.  All  of  these  mechanisms  require  additional  energy  input  into 
the  system  and  may  explain  the  apparent  toughening  of  the  material., 

Cohesive  zone  models  based  on  the  Dugdale-Barrenblatt  approach  of  placing 
a  closing  pressure  on  the  faces  of  a  discrete  crack  could  be  used  t-6  account  for 
the  nonlinear  effects  of  the  fracture  process  zone  (FPZ) .  Computer  codes  based 
on  such  concept  have  been  developed. is  generally  assumed  that  the  closing 
pressure  is  a  function  of  the  crack  opening  displacement  (w)  Since  the  results 
of  finite  element  computations  are  quite  sensitive  to  the^hape  of  the  assumed 
closing  pressure  vs.  w  relationship  it  is  essential  to  determine  such 
relationship  as  accurately  as  possible. 

Normally  closing  pressure  vs.  w  relationships  are  calibrated  by  comparing 
experimentally  measured  macroscopic  displacements  (such  as  load- joint  deflection 
in  beam  test)  with  those  calculated- numerically .  Such  a  procedure  does  not  give 
a  precise  or  an  unique  relationship  between  closing  pressure  vs.  crack  opening 
displacement. 

In  this  study,  crack  profiles  as  well  as  the  full  field  microscopic 
displacements  on  the  fracture  process  zone  were  accurately  measured  by  laser 
holographic  interferometry^  With  this  technique  it  was  possible  to  measure  the 
two-dimensional  displacement  field  around  a  crack  tip  with  an  accuracy  of  a 
quarter  of  a  micron.  "Such  accurate  measurement  of  strain  field  in  FPZ  provided 
a  be tter'"Ohder standing  of  the  fracture  process  in  concrete. 

Iphe  first  two  parts  of  this  report  discuss  mode  I  fracture  of  concrete 
which  was  studied  using  center  notched  plate  specimens.  Initially,  crack 
profiles  were  measured  by  illuminating  the  specimens  only  from  a  single 
direction.  This  was  sufficient  to  obtain  crack  opening  displacements.  To 
obtain  2-D  strain  field  three  simultaneous  holograms  were  made  by  illuminating 
the  specimens  from  j  different  directions.  To  analyze  three  different 
simultaneous  holographic  patterns,  a  digital  image  analysis  technique  was 
developed.  The  results  from  holographic  study  were  compared  with  finite  element 
analysis  based  on  linear  elastic  fracture  mechanics.  This  comparison  permitted 
a  new  definition  of  FPZ  as  well  as  more  accurate  relationship  between  closing 
pressure  and  crack  opening  displacement. 

-In  the  lasfpart  of  this  report,  a  study  of  mixed  mode  crack  propagation 
under  compression  is  described.  Mortar  plates  were  subjected  to  uniaxial 
compression.  These  specimens  had  notched  which  were  inclined  from  the  direction 
of  compression.  The  crack  trajectory  and  crack  opening  and  sliding  displacement 
were  measured  with  laser  holographic  interferometry.  The  results  were  compared 
with  LEFM  solutions  obtained  using  finite  element  analysis.  A  cohesive  crack 
model  which  includes  both  normal  and  shear  traction  was  developed  for  the  case 
of  mixed  mode  crack  propagation. 
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Part’  I. 


COHESIVE  CRACK  MODELS  EXAMINED  WITH  LASER  HOLOGRAPHIC  MEASUREMENTS 


1 . 1  Introduction 

Conventional  linear  elastic  fracture  mechanics  cannot  explain  fracture  in 
concrete  and  mortar.  Unlike  traditional  LEFM  materials  which  show  a  constant 
fracture  toughness,  concrete  and  mortar  toughen  (show  more  resistance  to 
fracture)  as  cracks  propagate  until  an  instability  results  and  the  structure 
fails.  This  increase  in  toughness  is  thought  to  occur  because  of  many  possible 
complex  mechanisms  which  take  place  during  fracture  in  concrete  and  mortar. 
These  mechanisms  include:  formation  of  a  zone  of  microcracking  around  the  crack 
tip;  crack  bridging,  where  aggregate  particles  span  the  crack  and  act  as 
ligaments  to  hold  the  crack  together;  and  crack  deflection,  where  the  crack 
changes  direction.  All  of  these  mechanisms  require  additional  energy  input  into 
the  system  and  may  explain  the  apparent  toughening  of  the  material. 

1.2  Cohesive  Crack  Models 

Cohesive  zone  models  based  on  the  Dugdale-Barrenblatt  approach  of  placing 
a  closing  pressure  on  the  faces  of  a  discrete  crack  (Fig.  I.lA)  could  be  used 
to  account  for  the  non-linear  effects  of  the  FPZ  and  aggregate  interlock. 
Hillerborg,  et.  al.[11  proposed  a  cohesive  model  for  cement  based  materials. 
This  model  assumes  that  the  closing  pressure  is  a  function  of  the  crack  opening 
(w) ;  Fig.  I.  IB).  Jenq  and  Shaht21  proposed  a  model  where  Kx  was  not  equal  to 
zero  (Fig.  I. 1C).  Cook,  Fairbanks,  Lawn  and  Mai131  also  used  a  model  where  K: 
was  not  equal  to  zero  and  which  had  a  degree  of  crack  length  dependence  (This 
model  was  for  ceramic  materials,  but  may  be  applicable  to  cement  based 
materials;  Fig.  I. ID).  It  is  assumed  in  the  Cook  model  that  the  material  has 
randomly  spaced  "ligaments"  which  cause  the  closing  pressure.  Short  cracks 
intersect  no  ligaments  so  there  Is  no  closing  pressure  while  long  cracks 
intersect  so  many  ligaments  that  a  steady  state  closing  pressure  is  used. 
Roelfstra  and  Wittmann141  and  Alvarado  and  Torrent151  argued  that  the  shape  of  the 
closing  pressure  vs.  w  curve  was  very  important  and  affected  the  predicted  post¬ 
peak  branch  of  the  load  vs.  crack  mouth  opening  displacement  (CMOD)  curve. 
Roelfstra  and  Wittmann  also  suggested  that  different  closing  pressure  vs.  w 
curves  were  needed  for  concrete  and  mortar  with  different  maximum  aggregate 
sizes.  Liaw,  Jeang,  Hawkins  and  Kobayashi161  proposed  closing  pressure  vs.  w 
relationships  based  on  concrete  strengths. 

1 . 3  Review  of  the  Experimental  Data 

The  crack  profile  data  used  during  the  course  of  this  research  was 
obtained  by  the  use  of  laser  holographic  interferometry.  A  brief  summary  of  the 
experimental  method  is  presented  here  and  further  detail  can  be  found  in 
Reference  10. 

A  center  notched  plate  specimen,  made  of  mortar,  was  tested  (Fig.  1.2). 
The  mortar  had  a  ratio  of  cement/sand/water  of  1/2.65/0.65  and  the  maximum 
particle  size  was  3mm  (0.125  inch).  The  specimens  were  loaded  in  parallel  with 
strain  gaged  steel  bars  to  provide  stable  crack  growth  (Fig.  I . 3A) . 

At  approximately  0.5  kN  (100  pound)  increments,  the  specimen  was 
hologramed  and  holograms  from  two  different  load  steps  were  combined  using  a 
technique  called  sandwich  holography171.  When  the  hologram  is  viewed,  the 
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surface  of  the  object  is  covered  with  fringes  and  the  crack  appears  as  a 
discontinuity  in  the  fringe  pattern  (Fig.  I.3B).  If  the  crack  is  crossed,  some 
fringe  orders  will  be  missing.  The  number  of  missing  fringe  orders  can  be  found 
by  tracing  a  line  around  the  crack  and  determining  the  number  of  fringes  crossed 
by  the  line.  The  opening  of  the  crack  is  proportional  to  the  fringe  counts. 

Miller,  Shah  and  Bjelkhagen181  showed  that  short  cracks  have  profiles  which 
are  basically  the  same  as  LEFM  profiles .  Longer  cracks  have  profiles  which  are 
much  thinner  than  LEFM  crack,  suggesting  that  applying  a  closing  pressure  to  the 
LEFM  profile  may  provide  the  correct  final  profile  (Fig.  1.4). 

1.4  Closing  Pressures  to  be  Studied 

Four  different  closing  pressures  have  been  chosen  for  study.  In  choosing 
these  closing  pressures,  consideration  was  given  to  both  the  similarities  and 
differences  between  the  experiments  which  produced  these  relationships  and  the 
holographic  experiment  which  will  be  used  here  for  comparison.  In  general, 
three  things  were  considered:  type  of  material  used  (i.e.  maximum  aggregate 
size,  w/c  ratio),  type  of  specimen  and  method  of  crack  width  measurement. 
Figure  1.5  is  a  summary  of  the  closing  pressure  relationships  studied  and  Table 
1.1  compares  available  properties. 

The  first  is  an  exponential  closing  pressure  proposed  by  Gopalaratanam  and 
Shah181.  This  closing  pressure  was  developed  from  tests  of  three  different 
materials;  a  concrete,  a  mortar  and  cement  paste. 

The  second  relationship  was  a  bilinear  curve  proposed  by  Roelfstra  and 
Wittmann151 .  As  was  previously  noted,  Roelfstra  and  Wittmann  thought  that 
maximum  aggregate  size  was  an  important  parameter  and  they  produced  a  series  of 
bilinear  curves  for  different  types  of  concrete  and  mortar.  The  bilinear  curve 
chosen  for  this  study  was  a  curve  for  a  mortar  which  had  the  same  maximum 
aggregate  size  as  the  mortar  used  in  the  holographic  study. 

A  trilinear  closing  pressure  proposed  by  Liaw,  Jeang,  Hawkins  and 
Kobayashi171 .  This  closing  pressure  was  chosen  for  analysis  because  the  crack 
line  wedge  loaded  double  cantilever  beam  specimen  used  by  Liaw,  et  al.  ,  was  very 
similar  to  half  the  center  cracked  plate  used  in  the  holographic  study  and 
because  the  closing  pressure  relationship  uses  strength  of  the  concrete  as  a 
parameter. 

The  final  closing  pressure  relationship  studied  was  a  linear  relationship 
proposed  by  Cedolin,  Iori,  and  Dei  Polit91.  The  importance  of  the  linear  closing 
pressure  (and  the  reason  for  its  inclusion  in  this  study)  is  that  this  is  the 
only  other  published  study  the  authors  could  find  which  used  an  optical  laser 
method  to  measure  crack  opening. 

1.5  Methodology  of  Analysis 

The  analysis  in  this  study  was  performed  by  using  a  finite  element  program 
where  the  given  closing  pressure  was  applied  to  the  crack  faces.  The  resulting 
crack  profile  was  then  compared  with  holographically  measured  crack  profiles . 
An  iterative  procedure  is  needed  because  the  closing  pressure  depends  on  the 
unknown  crack  opening,  w.  If  the  closing  pressure  is  too  high,  the  calculated 
crack  opening  displacements  become  negative  (overlap  of  crack  faces) .  Special 
provisions  have  to  be  made  in  the  finite  element  program  to  prevent  overlapping. 

1.6  Material  Properties 

In  addition  to  the  plate  specimens,  50.8  mm  X  101.6  mm  (2  in  X  4  in) 
cylinder  specimens  of  the  same  mortar  were  tested  in  uniaxial  compression  at 
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comparable  age.  Two  linear  variable  differential  transducers  (LVDTs)  with  a 
50.8  mm  (2  in)  gage  length  were  mounted  on  diametrically  opposite  sides  to 
measure  axial  deformation.  The  initial  tangent  modulus  was  used  during  the 
analysis . 

Measured  tensile  strength  shows  considerable  scatter  and  there  is  a  marked 
variation  depending  on  the  method  used  to  determine  tensile  strength.  Depending 
on  the  method  used,  tensile  strength  may  vary  from  7  to  14%  of  the  compressive 
strength  (f'c)t101.  This  scatter  will  affect  the  consistency  of  the  model.  In 
this  study,  tensile  strength  was  assumed  to  be: 

f't  -  0 . 10f'e 

1.7  Analysis  of  the  Closing  Pressure  Models  for  Kx  -  0 

The  applied  load  Pa  and  Kx  -  0  were  specified  and  then,  the  crack  length 
was  calculated.  This  was  done  by  applying  the  experimental  load  to  the  mesh 
and  then  propagating  the  crack  3  mm  (1/8  in)  at  a  time.  At  each  increment  of 
crack  length  the  closing  pressure  was  applied  and  a  value  of  Kx  calculated  until 
a  value  of  Kj  —  0  was  found. 

For  some  closing  pressures  it  was  not  possible  to  find  a  Kx  -  0  crack 
length  since  became  either  constant  or  began  to  increase  (Fig.  1.6).  This 
occurred  only  at  the  peak  load.  In  this  case  the  minimum  Kx  was  used  for  the  Kx 
-  0  length,  but  these  cases  will  be  noted  in  the  analysis. 

Figure  1.7  shows  the  predicted  crack  length  for  Kj-0  vs.  measured  crack 
length.  Except  for  the  trilinear  closing  pressure,  all  the  closing  pressures 
predict  the  correct  crack  length  within  the  experimental  accuracy  (shown  by  the 
dashed  lines).  Therefore,  accurate  prediction  of  thee  crack  length  is  not 
restricted  to  a  unique  closing  pressure. 

A  better  evaluation  of  the  different  closing  pressures  results  from 
comparing  the  experimental  crack  profiles  vs.  predicted  crack  profiles.  In 
order  to  provide  an  objective  comparison,  an  error  definition  is  proposed: 

s  (wko  *  wh> 
error  -  - 

Swh 


where:  wh  -  measured  holographic  crack  opening  at  a  node 

wK0  -  crack  opening  for  Kt  -  0  profile  at  the  same  node 

The  sum  is  taken  over  the  number  of  finite  element  nodes  along  the  crack  face 
(spaced  at  3mm  apart).  A  positive  error  means  the  closing  pressure  was  not 
large  enough  to  close  the  LEFM  crack  profile  down  to  the  holographic  profile 
while  a  negative  error  means  the  closing  pressure  was  too  strong  and  closed  the 
LEFM  crack  profile  down  too  much. 

Figure  1.8  shows  the  error  for  all  the  specimens  tested.  From  figure  1.8 
it  is  apparent  that  comparable  results  were  obtained  for  the  cracks  developed 
on  specimens  of  same  strength  (i.e.  same  results  for  specimens  tested  on  1/30 
and  1/31)  which  may  indicate  a  strength  dependance  on  the  closing  pressure. 
Also,  note  that  with  the  bilinear  closing  pressure,  short  cracks  are  closed  too 
much  even  though  for  the  longer  cracks  there  is  good  agreement  between 
experimental  and  predicted  results.  This  may  indicate  the  use  of  less  closing 
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pressure  near  the  crack  tip  (most  important  for  short  cracks)  while  for  longer 
cracks  this  closing  pressure  may  be  adequate.  This  is  not  surprising  since  it 
was  previously  shown  that  short  cracks  are  nearly  LEFM  cracks  and  no  closing 
pressure  is  required.  In  the  other  hand,  the  trilinear  closing  pressure  shows 
a  somewhat  opposite  behavior.  For  short  cracks  there  is  a  relatively  good  match 
while  for  long  cracks  the  closing  pressure  appears  to  be  too  strong. 

1.8  Proposed  Adjustments  to  the  Bilinear  Closing  Pressure 

Out  of  the  four  closing  pressures  examined,  the  bilinear  closing  pressure 
provided  the  "best"  fit,  although  the  fit  was  not  very  good.  The  next  step 
seemed  to  be  to  see  if  the  parameters  of  the  bilinear  curve  could  be  logically 
adjusted  to  provide  a  better  fit.  Since  the  crack  lengths  and  holographically 
measured  crack  openings  are  known,  it  is  possible  to  calculate  the  required 
closing  pressure  distribution  that  exactly  matches  the  experimental  results. 
Figure  1.9  shows  some  sample  calculated  closing  pressure  distributions. 

Cook,  et.  al.  [4]  proposed  that  crack  length  dependence  could  be  accounted 
for  in  the  cohesive  crack  model  by  assuming  that  there  is  some  initial  distance 
over  which  the  closing  pressure  does  not  act  (see  introduction).  This  agrees 
with  figure  1.4  that  shows  good  agreement  between  experimental  and  LEFM  crack 
profiles  (no  closing  pressure)  for  short  cracks.  This  model  can  be  incorporated 
into  this  study  by  releasing  the  closing  pressure  from  some  of  the  nodes,  i.e. 
making  some  nodes  ahead  of  the  original  not  tip  traction  free. 

An  objective  procedure  to  obtain  the  best  bilinear  closing  pressure  for 
all  the  cracks  analyzed  was  proposed.  The  penalty  function  f  given  by: 

f  -  S  errort2 


where : 

errort  -  relative  error  for  load  step  i 

was  defined  and  a  bilinear  closing  pressure  defined  by  four  parameters  (ft, 
flt  wx  and  w0;  Fig.  1.9  )  is  considered.  An  algorithm  was  applied  to  minimize 
the  penalty  function  for  each  specimen.  The  same  algorithm  was  applied  for 
various  lengths  of  the  traction  free  zone  in  front  of  the  notch  tip  (as  proposed 
be  Cook)  with  the  best  results  for  a  traction  free  zone  of  8  mm  (1/8  in)  .  A 
proposed  bilinear  closing  pressure,  based  on  the  optimization  procedure  is  shown 
in  figure  1.9. 

Table  1 . 2  shows  the  results  of  the  optimization  precess .  Even  though  only 
two  different  strengths  were  tested  some  conclusions  can  be  made.  There  is  a 
clear  strength  dependence  on  the  stress  parameters  (ft  and  fx)  while  the  crack 
opening  parameters  converged  to  the  same  solution  regardless  of  the  strength  of 
the  material.  Also,  note  that  for  the  specimens  tested,  the  optimized  value  of 
ft  is  about  10.5%  of  f'e  which  agrees  with  previous  experimental  measurements  of 
the  tensile  strength  of  concrete. 

Figure  1.10  shows  the  relative  error  plots  for  the  proposed  bilinear 
closing  pressure  and  typical  crack  profile  comparison.  Note  that  there  is  a 
good  prediction  of  the  crack  opening  displacements  although  the  crack  lengths 
are  over  estimated  including  a  segment  close  to  the  crack  tip  with  either  small 
or  null  predicted  crack  opening. 

Jenq  and  Shah  suggested  that  K:  did  not  need  to  be  equal  to  0.  If  the 
closing  pressure  is  applied  Kx  corresponding  to  the  experimentally  observed  crack 
length  is  calculated,  an  upper  bound  for  each  specimen  is  observed  (Fig.  I. 11). 
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An  upper  bound  is  used  because  the  holograms  are  made  at  discrete  points  during 
the  loading  and  it  is  not  possible  to  tell  if  the  crack  is  just  about  to 
propagate,  has  just  propagated  or  is  somewhere  in  between.  Figure  1.12  shows 
the  relative  error  plots  and  typical  crack  profiles  corresponding  to  the  upper 
bound  value  of  Kx  from  figure  1. 11.  Note  that  there  is  no  significant  effect  on 
the  results. 

1.9  Discussion  and  Conclusions 

This  study  clearly  shows  that  the  cohesive  crack  model  can  correctly 
predict  crack  behavior  but  that  the  assumption  that  the  closing  pressure  depends 
only  on  the  crack  opening,  w,  may  not  be  correct.  The  results  show  that 
strength,  maximum  particle  size,  and  crack  length  are  also  important  parameters. 
The  assumption  that  Kt  -  0  seems  to  be  acceptable  as  long  as  the  actual  Kx  -  0 
crack  length  is  calculated  and  not  assumed.  In  addition,  it  is  very  important 
to  note  that  the  values  of  modulus  and  tensile  strength  greatly  affect  the 
model.  Due  to  the  limited  amount  of  accurate  data  on  crack  lengths  and 
crack  opening  displacements,  the  authors  cannot  yet  propose  exact  modifications 
to  the  closing  pressure  functions,  although  the  Cook,  et.  al.t4J  model  seems  to 
be  a  good  start  on  including  crack  length  effects  into  the  cohesive  crack  model. 
The  authors  hope  that  as  more  holographic  experiments  are  performed  on  different 
mixes  and  specimens,  a  clearer  picture  of  the  necessary  changes  to  the  cohesive 
zone  model  may  emerge. 

This  research  did  not  address  the  question  of  making  the  area  under  the 
closing  pressure  curve  equal  the  fracture  energy,  Gf.  While  the  fracture  energy 
is  an  important  parameter,  it  seems  that  there  is  no  agreed  upon  value  or 
function  for  Gt.  Values  of  G£  can  be  found  in  literature  to  justify  almost  any 
closing  pressure.  Therefore-,  more  work  needs  to  be  done  on  defining  Gf  for 
concrete  before  it  will  be  a  useful  parameter  for  closing  pressures. 
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TABLE  1.1 


COMPARISON  OF  CONCRETE  PROPERTIES 


REFERENCE/CLOSING 

MAXIMUM 

f'c 

f't 

PRESSURE 

AGGREGATE 

(mm) 

MPa 

MPa 

w/c 

Cedolin  (Linear) 

8 

41.4 

3.2 

.5 

Roelfstra( Bilinear) 

3 

N/A 

4.7 

.58 

Jeang  (Trilinear)* 

6-12 

27-62 

2.3-3.5f 

N/A 

Gopalaratanara  (Exp)* 

0-10 

33-47 

N/A 

.  4  - 

Miller  (Holographic) 

3 

28  &  47 

N/A 

.6 

Multiple  mixes  used  to  find  closing  pressure 
+  From  formula  f't  -  c  fit  to  experimental  data 

TABLE  1.2 


OPTIMIZED 

BILINEAR 

CLOSING  PRESSURE 

TEST 

f  * 

ft 

fi 

wi 

wo 

DATE 

(MPa) 

(MPa) 

(MPa) 

(micron) 

(micron) 

8 

28 

3.13 

1.45 

3.81 

73.66 

30 

47 

4.47 

2.16 

3.81 

67.56 

31 

47 

4.91 

1.75 

3.82 

69.60 
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Part  II.  MEASUREMENT  OF  STRAIN  FIELD  IN  FRACTURE  PROCESS  ZONE  OF  MORTAR 

11 . 1  Introduction 

The  study  presented  in  part  I  shows  that  the  cohesive  crack  model  approach 
can  be  used  to  predict  both  crack  extension  and  crack  opening  displacements. 
However,  the  fracture  process  in  concrete  is  not  limited  to  discrete 
macrocracks.  Several  phenomenon  have  been  associated  to  the  so  called  fracture 
process  zone  (FPZ)  such  as  microcracking  around  the  crack  tip  and  strain 
softening.  A  key  to  a  better  understanding  of  fracture  of  concrete  is  the 
accurate  observation  of  the  fracture  process  zone. 

The  first  observation  of  this  kind  was  made  by  Cedolin,  Dei  Poli  and  Iori. 
They  used  laser  moire  interferometry  for  wide  field  measurements  in  specimen 
under  tension.  White  light  moird  interferometry  was  also  used  by  Du,  Kobayashi 
and  Hawkins1111  to  measure  crack  opening  displacement  along  the  fracture  process 
zone  (FPZ)  for  concrete  specimens  under  three  point  bending.  Regnault  and 
Bruhwiler1121  used  moire  holographic  interferometry  to  measure  strain  field  on 
wedge  splitting  test  specimens  and  concluded  that  the  model  based  on  the  exact 
knowledge  of  traction  free  crack  tip  location  does  not  describe  the  material 
behavior  in  a  realistic  way.  Raiss,  Dougill  and  Newman1131  also  used  laser  moire 
interferometry  to  measure  strain  field  on  direct  tension  specimens. 

In  the  technique  described  in  this  report,  three  simultaneous  holograms 
were  taken  by  illuminating  the  specimen  from  three  different  directions.  In 
addition,  digital  image  analysis  techniques  were  developed  to  accurately 
evaluate  wide  field  in  plane  deformation.  The  results  of  the  strain  field 
calculation  obtained  using  this  newly  developed  technique  for  mortar  center 
notched  plate  specimens  are  discussed. 

11. 2  Experimental  Program 

Center  notched  mortar  plate  specimens  were  used  (Fig.  II. 1).  3  in  x  6  in 
(  76  mm  x  152  mm)  cylinders  were  also  cast  in  order  to  determine  the  compressive 
strength  and  initial  modulus  of  elasticity.  Two  cylinders  with  1/2  in  (13  mm) 
strain  gages  were  tested  in  a  closed  loop  testing  system  under  strain  control 
the  day  after  the  center  notched  specimen  test  was  performed. 

Holographic  interferometry  was  used  to  measure  crack  opening  displacements 
and  the  strain  field  on  a  2  in  x  3  in  (51  mm  x  76  mm)  area  in  front  of  each 
notch  (i.e.,  top  and  bottom;  Fig.  II. 1). 

11. 3  Holographic  Interferometry  Test 

Each  interference  pattern  represents  the  component  of  the  displacements 
onto  the  sensitivity  vector  kt71 .  In  the  optical  set-up  shown  in  figure  II. 2 
all  the  optical  paths  are  on  a  plane  of  constant  height  (y-constant) . 
Therefore,  the  sensitivity  of  the  optical  system  to  vertical  displacements  is 
small.  In  addition,  the  loading  configuration  restrains  displacements  in  the 
y  direction.  Thus,  two  displacement  components  are  sufficient  to  describe  to 
motion  of  any  point  on  the  specimen,  namely  the  in  plane  component  (x)  and  the 
out  of  plane  component  (z) .  In  order  to  obtain  a  two  dimensional  displacement 
field  (in  plane  ana  out  of  plane),  in  principle,  two  interferograms  are 
required.  However,  since  there  is  no  sign  associated  to  the  fringe  counts  a 
third  interferogram  is  necessary.  Figure  II.  3  is  an  schematic  view  of  the 
optical  set-up  with  three  illumination  directions  (R,L  and  C)  used  for  this 
s  tudy . 
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For  each  load  step  a  different  sandwich  hologram1141  is  recorded  from  each 
illumination  direction  (holograms  R,L  and  C) .  Due  to  the  incremental  nature  of 
the  holographic  measurements,  it  is  necessary  to  make  a  series  of  sandwich 
holograms  at  continuous  load  levels.  If  a  load  increment  is  not  recorded  none 
of  the  future  deformation  states  can  be  computed. 

11. 4  Enhancement  of  Holographic  Interferograms  by  Image  Analysis 

The  holographic  interferograms  were  acquired  into  a  Tracor  Northern  image 
analysis  system  using  a  slow  scan  high  resolution  COHU  camera.  A  schematic 
arrangement  for  reconstruction  of  the  holographic  images  and  acquisition  into 
the  image  analysis  system  is  show  in  figure  II. 3.  The  acquired  image  (grey 
image)  contains  a  range  of  intensities  (0-255)  representing  different  grey 
levels . 

The  bright  and  dark  bands  of  the  holographic  fringes  can  be  represented 
as  a  binary  image.  In  a  binary  image  the  pixel  intensities  assume  only  two 
values,  ON  (or  1)  for  perfect  white  and  OFF  (or  0)  for  perfect  black.  Due  to 
variations  of  intensity  of  the  object  beam  throughout  the  surface  of  the  object, 
it  is  possible  that  a  bright  band  at  a  section  of  the  image  has  a  lower 
intensity  magnitude  than  the  dark  bands  at  other  segments  of  the  image. 
Standard  methods  of  image  processing  such  as  thresholding  from  a  certain  pixel 
intensity  are  thus  inapplicable.  This  situation  is  represented  by  the  plot  in 
figure  II. 4a  which  shows  the  intensity  variations  along  a  cross  section  of  a 
holographic  grey  image. 

If  Ix  and  I2  are  the  intensities  of  the  individual  images  in  a  sandwich 
hologram,  the  intensity  of  the  resulting  interference  pattern  is: 

I  -  (Ix  +  I2)  +  2(IX  la)1'2  cos  A*  (1) 

From  this  equation  it  is  apparent  that  the  intensity  value  of  any  pixel  can  be 
separated  into  two  terms.  One  is  the  average  intensity  of  the  image  (Ix  +  I2)  , 
and  the  second  depends  on  the  phase  difference  between  the  two  wave  fronts  (A^)  . 
The  separation  of  this  two  terms  can  be  achieved  through  the  use  of  the  image 
analysis  system. 

Assume  that  the  grey  image  of  the  sandwich  holography  interferogram  is 
represented  as  IMAGE_1  where  as  IMAGE_2  is  the  average  of  the  independent  wave 
fronts  acquired  separately  (Fig.  II. 4a).  Subtraction  of  IMAGE_2  from  IMAGE_1 
results  in  the  isolation  of  the  interferometric  effect,  represented  as  IMAGE_3 
in  figure  II. 4b.  Since  negative  values  are  not  allowed  as  meaningful  intensity 
values,  IMAGE_3  results  in  zero  for  all  the  pixels  on  dark  bands.  The  bright 
bands  retain  a  positive  magnitude.  By  selecting  all  positive  values  a  binary 
image  is  created  (Fig.  II. 4c). 

11. 5  Fringe  Count 

IMAGE_3  is  the  enhanced  interferometric  image  in  which  the  intensity  value 
is  zero  for  dark  bands  and  one  for  bright  bands.  Further  image  processing 
allows  to  assign  intensity  values  to  each  fringe  corresponding  to  the  same 
fringe  order. 

Figure  II. 5  shows  a  fully  processed  interferogram  in  which  same  intensity 
values  correspond  to  same  fringe  orders.  Note  that  there  is  a  different 
intensity  level  for  each  half  a  fringe  (i.e.,  different  intensity  for  dark  and 
bright  bands).  The  fringe  count  between  any  two  points  can  now  be  easily 


obtained  by  computing  half  the  difference  of  the  corresponding  intensity  levels 
(e.g.,  fringe  count  between  A  and  B  in  figure  II. 5  is  (90  -  65)/2  -  12.5). 

Even  though  the  use  of  the  image  analysis  system  results  in  very  fast 
fringe  counts,  speed  is  not  the  main  advantage  of  this  technique.  The  isolation 
of  the  interferometric  effect  on  the  holographic  images  and  the  binarization 
process  result  in  accurate  and  objective  fringe  counts.  The  gradual  change  in 
intensity  from  a  dark  to  a  bright  band  make  it  very  difficult  to  obtain 
consistent  fringe  counts  when  conventional  (by  hand)  methods  are  used. 

11. 6  Evaluation  of  Holograms 

The  holodiagramtl5J  shown  in  figure  II. 6  is  a  simple  geometrical  tool  for 
computation  of  displacements  based  on  interferometric  readings.  In  figure  II. 6 
the  distance  I-P-0  from  Illumination  point  to  object  to  observation  point 
(usually  taken  as  the  position  of  the  holographic  plate)  remains  constant  if  the 
object  (P)  moves  along  ellipse  Ex  (or  an  ellipsoid  in  the  most  general  three 
dimensional  case)  with  foci  at  illumination  and  observation  points.  Now,  define 
ellipse  E2  such  that  the  distance  I-P-0  increases  by  A/2  (180°  phase  difference) . 
An  interferogram  of  the  object  that  moved  from  any  point  on  ellipse  E2  to  any 
point  on  ellipse  E2  will  show  a  dark  band  (half  a  fringe)  due  to  destructive 
interference  of  the  wave  fronts.  In  general,  the  number  of  fringes  is 
proportional  to  the  increase  (or  decrease)  of  the  distance  I-P-0. 

From  figure  II.  6  it  is  apparent  that  there  are  an  infinite  number  of 
possible  displacement  vectors  which  result  on  the  same  fringe  count.  The  number 
of  possible  displacement  vectors  is  reduced  to  four  by  using  two  holograms  with 
different  illumination  points  and  fringe  counts  nx  and  n2  (inset  of  figure  II. 6). 
Using  a  fringe  count  from  a  third  hologram,  it  can  be  determined  if  nx  and  n2 
have  the  same  or  opposite  sign. 

11. 7  Crack  Profiles. 

For  the  computation  of  the  crack  opening  displacements  (COD)  the  modified 
Nelson  and  McCrickerd  1161181  method  was  used.  This  method  assumes  the  out  of 
plane  motion  to  be  negligible  (i.e.  COD  is  the  main  component  of  the 
displacement)  and  that  displacements  are  small.  Since  the  direction  of  motion 
is  assumed,  one  hologram  is  sufficient.  This  was  described  in  part  I. 

An  alternative  solution  is  to  assume  a  two  dimensional  displacement 
solution.  The  displacement  is  computed  as  the  difference  in  position  of  the 
intersection  of  the  ellipses  from  two  illumination  directions.  The 
holographic  fringe  counts  can  be  combined  in  three  possible  ways  (  R  and  L,  R 
and  C,  C  and  L) .  As  mentioned  before,  the  third  fringe  count  is  necessary  to 
select  the  correct  value  among  the  four  possible  displacement  solutions  for  each 
combination.  Figure  II. 7  shows  a  crack  profile  computed  using  all  three 
combinations.  Note  that  there  is  an  improvement  in  the  scatter  on  the  COD 
calculations  with  respect  to  the  Nelson  and  McCrickerd  results.  Also,  note  that 
taking  into  account  the  out  of  plane  contribution  to  the  fringe  counts  results 
in  smaller  CODs  even  for  an  experimental  set-up  designed  for  mode  I  crack 
propagation.  This  means  that  there  exist  a  measurable  component  of  the  crack 
face  displacement  in  the  out  of  plane  direction  (mode  II)  which  can  explain  the 
presence  of  compressive  strain  behind  the  crack  tip  as  reported  by  Bruhwiler. 
The  doted  line  in  figure  II. 7  represents  the  crack  sliding  displacement  (CSD) 
in  the  XZ  plane. 

11. 8  Strain  Fields. 
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Fringe  counts  relative  to  the  left  face  of  the  crack,  or  the  symmetry  line 
for  point  beyond  the  crack  tip,  were  taken  at  every  point  on  a  3  in  x  3  in  (76 
mm  x  76  mm)  grid  every  1/8  in  (3  mm)  for  all  three  holographic  images  (R,L  and 
C)  .  A  computer  program  was  used  to  calculate  the  relative  displacement  of  every 
point  relative  to  the  left  face  of  the  crack  by  locating  the  intersection  point 
of  the  holodiagrams  corresponding  to  each  illumination  direction.  The  strain 
field  was  obtained  by  differentiation.  Figures  II.  8  show  the  strain  field 
(e**)  around  the  bottom  crack  for  crack  lengths  of  1.375  in  (35  mm),  1.875  in  (48 
mm)  and  2.25  in  (57  mm),  corresponding  to  1324  lb  (6.0  KN) ,  1461  lb  (6.6  KN)  and 
1601  lb  (7.2  KN)  of  applied  load.  It  is  important  to  notice  that  during  the 
experimental  program,  before  any  holograms  were  made,  the  specimens  were  placed 
under  a  800  lb  (3.6  KN)  preload.  This  preload  is  necessary  to  keep  the  specimen 
optically  stable.  The  strain  fields  shown  in  figures  II. 8  include  only  the 
measured  additional  strain  after  the  initial  preload  state. 

In  general,  a  high  strain  (>100 /ie)  region  can  be  found  around  and  in  front 
of  the  crack  tip  with  a  maximum  value  of  approximately  300  ^e.  Behind  the  crack 
tip  there  is  a  zone  of  negative  strain  which  suggests  relaxation  of  tensile 
strain  with  respect  to  the  initial  preload  stage. 

II.  9  Definition  of  the  Fracture  Process  Zone 

The  zone  where  the  measured  strain  values  exceed  100  fie  (limit  to  define 
nonlinearity)  are  shown  in  figures  II. 9  for  three  crack  lengths.  Also  shown  in 
the  figure  are  the  zones  behind  the  crack  tip  where  tensile  strain  was  detected. 
Note  that  the  end  of  these  zones  correspond  approximately  to  the  tensile  strain 
value  of  zero.  There  appears  to  be  a  traction  free  zone  of  approximately  0.25 
in  front  of  the  original  notch  tip  regardless  of  the  crack  length  and  COD 
indicated  by  wc  in  figures  II. 9).  This  agrees  with  the  model  suggested  by  Cook, 
et.  al.t3J 

An  alternate  approach  to  define  the  FPZ  is  to  find  the  deviations  from  the 
linear  elastic  solution.  The  specimen  was  analyzed  using  a  finite  element 
program.  In  order  to  make  a  valid  comparison  with  the  experimental  results, 
figures  II. 10  show  the  LEFM  strain  distribution  corresponding  to  the  difference 
between  the  given  load  level  and  the  initial  preload.  Note  that  the  strain 
levels  surrounding  the  crack  tip  exceed  those  from  the  holographic  measurements 
(Fig.  II. 8) 

Behind  the  crack  tip  the  LEFM  solution  shows  sharp  unloading.  At  1601  lb 
the  strain  drops  from  500  pe  at  the  crack  tip  to  zero  within  1/8  in  behind  the 
crack  tip.  In  contrast,  the  observed  strain  relaxation  shown  on  figure  II. 8 
follows  a  gradual  descent. 

In  figures  II. 11,  the  differences  in  strain  fields  between  the  LEFM 
solution  and  the  experimentally  measured  values  are  shown.  It  is  arbitrarily 
assumed  that  the  differences  are  significant  when  they  exceed  the  strain  value 
of  60  /i€ .  The  zones  marked  A  in  figure  II. 11  are  the  regions  where  the  LEFM 
solution  predicts  tensile  strain  values  60  fit  higher  than  those  measured 
holographically.  Note  that  this  zone  remains  essentially  constant  regardless 
of  crack  length.  Since  this  zone  is  relatively  small  and  does  not  change  with 
crack  length  a  model  based  on  modified  LEFM  may  be  a  possible  approach.  This 
is  discussed  in  the  next  section.  The  zones  B,  with  negative  values,  typically 
found  behind  the  crack  tip  show  that  the  relaxation  of  the  observed  strain  is 
60  pe  less  than  that  in  the  elastic  solution.  Zones  B  can  be  defined  as  the 
wake  of  the  fracture  process  zone  (WFPZ)  and  they  enclose  the  area  where 
extensive  microcracking  has  been  developed  and  tensile  forces  are  still 
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transmitted  through  the  crack.  Behind  this  zone,  the  fracture  process 
approaches  a  traction  free  condition. 

11. 10  Closing  Pressure  Analysis 

The  cohesive  type  of  model  with  a  bilinear  closing  pressure  described  in 
part  I  was  used  to  simulate  the  effect  of  the  fracture  process  zone.  To  examine 
whether  such  a  model  can  accurately  predict  the  measured  strain  field,  the  same 
finite  element  program  with  singular  elements  was  used.  The  predicted  crack 
profiles  in  figure  11.12  show  a  good  correlation  with  the  experimental  values. 

Figures  11.13  show  the  difference  in  strain  field  predicted  using  the 
closing  pressure  approach  and  the  observed  values.  Note  that  there  is  a  good 
correlation  of  the  results  and  the  zones  with  significant  deviation  from  the 
experimental  results  are  practically  eliminated  (±  60  fie). 

11. 11  Conclusions 

a.  Sandwich  hologram  interferometry  with  multiple  sensitivity  vectors  can  be 
used  to  measure  crack  opening  displacements  and  strain  fields  in  mortar 
specimens  under  tensile  stress.  When  holographic  interferometry  is  used 
to  measure  Mode  I  crack  opening  displacements  the  effect  of  out  plane 
motion  should  be  taken  into  account. 

b.  Digital  image  analysis  facility  provided  a  faster,  more  consistent  and 
more  accurate  method  of  fringe  count  than  the  manual  method. 

c.  There  is  a  region  of  high  tensile  strain  that  moves  with  the  tip  of  the 
propagating  crack.  The  material  behind  the  crack  tip  experiences  strain 
relaxation.  However,  the  gradual  nature  of  the  strain  relaxation 
demonstrates  the  existence  of  tensile  force  transmitted  through  the  crack 
faces. 

d.  Zones  of  nonlinear  behavior  can  be  located  by  computing  strain  field 
deviations  form  the  linear  elastic  solution. 

e.  A  bilinear  closing  pressure  vs.  COD  cohesive  crack  model  has  been  proposed 
to  predict  both  crack  profiles  and  strain  field,  fully  characterizing  the 
material  behavior. 
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Part  III.  A  STUDY  OF  MIXED  MODE  CRACK  PROPAGATION  IN  MORTAR  USING  HOLOGRAPHIC 
INTERFEROMETRY 

111 . 1  Introduction 

In  a  mode  I  pure  tensile  situation  described  in  parts  I  and  II,  the 
application  of  a  cohesive  crack  model  is  beginning  to  be  understood1171 1185 .  Most 
practical  situations  however  involve  a  combined  presence  of  tensile  and  shear 
cracking.  This  creates  a  mixed  mode  condition  involving  both  mode  I  and  mode 
III  stress  intensity  factors.  It  is  important  to  study  if  the  theoretical 
developments  of  mode  I  condition  can  be  extended  in  principle  to  a  mixed  mode 
situation.  Only  a  limited  number  of  studies  have  been  conducted  on  mixed  mode 
loading  of  concrete  and  other  quasi-brittle  materials  [i9H20]  [21]  [22]  [23]  [24] .  ^ost 
of  this  work  is  under  tensile  shear  condition.  There  is  a  lack  of  information 
about  compression  shear  crack  initiation  which  is  being  addressed  in  this 
report.  It  is  also  important  to  study  the  phenomena  involved  when  the  relative 
magnitude  of  the  mode  I  and  mode  III  stress  intensity  values  are  varied  at  the 
crack  tip  creating  different  types  of  mixed  mode  situations.  This  is 
particularly  important  because  of  the  coupling  between  the  shear  and  tensile 
openings  of  crack  in  such  materials  where  the  aggregate  interlocking  and  crack 
surface  irregularities  could  affect  the  final  crack  profile. 

111. 2  Theoretical  Background 

The  theoretical  considerations  are  based  on  an  infinite  plate  under 
uniaxial  compressive  stress  in  the  far  field  (Fig.  III.l).  The  plate  has  an 
existing  crack  length  at  an  angle  /9  to  the  direction  of  the  applied  load.  The 
crack  can  be  considered  to  be  subjected  to  the  far  field  stress  condition  shown 
in  the  figure  III.l  where: 

oy  -  -a  sin0,  ax  -  -a  cos/9,  -  -a  sin0  cos/9 

where  for  compressive  stress  a  >  0  the  mode  I  and  mode  III  stress  intensity 
factors  at  the  crack- tip  are  given  by: 

Kj  -  -a  7  ira  sin2/9,  KtI  -  -a  J  1 ra  sin/0  cos/3  (o>0) 

It  is  hence  possible  to  vary  the  ratio  of  Kj/Kn  by  varying  the  inclination 
angle  /9.  Negative  value  of  Kx  implies  closure  of  the  original  notch.  Negative 
value  of  KXI  does  not  have  any  such  significance  but  affects  the  direction  of 
crack  propagation.  Crack  initiation  of  the  inclined  crack  under  compression  is 
associated  with  closure  of  the  original  notch.  Cracks  are  initiated  under  a 
combination  of  compression  and  shear.  Note  that  once  the  crack  initiates  and 
propagates,  the  propagating  crack  is  under  tension  and  shear. 

The  equation  above  provide  an  insight  into  the  nature  of  the  problem  but 
are  valid  only  for  an  infinite  plate.  Due  to  the  finite  size  of  the  specimen  and 
the  gradual  curving  of  the  propagating  crack,  is  was  essential  to  conduct  an  FEM 
analysis  of  this  problem  at  different  stages  of  crack  propagation. 

111. 3  Finite  Element  Analysis 

In  order  to  study  the  Stress  Intensity  Factors  (SIFs)  and  displacements 
of  the  actually  tested  specimens  a  FEM  program  with  quarter  point  singular 
elements  was  usedC251.  Figure  III. 2  shows  a  typical  FEM  mesh  used  during  this 
study.  Note  the  triangular  elements  (singular  elements)  at  the  notch  tips.  As 
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the  crack  propagates  it  was  necessary  to  modify  the  FEM  mesh  accordingly.  This 
was  done  by  using  an  interactive  code  FRANC  which  has  automatic  remeshing 
capabilities.  In  the  experimental  program  four  types  of  specimens  were  testes 
which  were  of  the  same  dimensions  except  for  the  initial  notch  inclination 
angles  y9.  Four  different  values  of  /9  were  used  (18°,  36°,  54°  and  12°).  A 
program  was  written  to  generate  the  four  types  of  input  meshes  corresponding  to 
the  four  inclination  angles . 

Figure  III. 3  compares  the  solutions  obtained  from  the  FEM  analysis  to  the 
solution  from  the  infinite  plate  analysis  discussed  above  and  the  solution  for 
finite  geometry  obtained  by  Murakami1261.  The  plotted  data  are  values  of  Kx  at 
the  tip  of  the  original  notch  of  inclination  /3  when  the  applied  compressive 
stress  is  equal  to  the  stress  at  which  the  crack  starts  to  propagate.  The 
figure  shows  that  the  values  obtained  from  the  FEM  solution  are  very  close  to 
Murakami's  solutions.  The  infinite  plate  solutions  are  found  to  give  somewhat 
lower  results. 

III. 4  Experimental  Program 

Specimens  were  made  with  3: 1:0. 6  mixture  of  sand,  type  I  cement  and  water. 
The  central  portion  of  the  specimen  was  8"  X  6”  X  3"  (20.3  X  15.2  X  7.6  cm). 
Two  separate  end  blocks  were  cast  for  each  end  and  measured  1.5"  X  6"  X  3"  (3.8 
X  15.2  X  7.6  cm;  Fig.  III.l).  The  end  blocks  were  made  to  reduce  the  effect  of 
the  end  friction  from  the  loading  platen.  It  has  been  observed  that  this  simple 
method  is  the  best  way  to  ensure  an  uniform  uniaxial  compression  in  the  central 
portion  of  the  specimen1271.  The  specimens  were  cast  in  plexiglass  molds.  The 
end  pieces  were  cast  simultaneously  with  the  center  piece  by  placing  1/16"  (0.16 
cm)  plexiglass  separators  in  the  molds.  The  notch  was  2"  X  0.078"  (5.08  X  0.2 
cm)  and  was  cast  by  placing  a  steel  blade  in  the  mold  while  casting.  The  blade 
was  attached  to  an  aluminum  frame  which  was  fixed  to  the  side  of  the  mold.  The 
blade  was  attached  by  a  screw  and  a  nut  arrangement  which  made  it  possible  to 
change  the  angle  of  the  plate  and  cast  all  four  types  of  specimens  with  the  same 
set  of  molds.  The  molds  were  oiled  and  the  plates  were  greased  prior  to  pouring 
the  mixture.  The  blade  were  removed  four  hours  after  casting  when  the  mortar 
had  set  not  yet  harden.  Once  the  mortar  hardens  it  is  difficult  to  remove  the 
blades  without  damaging  the  specimen.  The  specimens  were  demolded  24  hours 
after  casting  and  cured  in  water  at  room  temperature  thereafter.  They  were 
tested  about  60  days  after  casting. 

The  specimens  were  tested  in  compression  under  120  Kip  MTS  loading 
machine.  They  were  capped  with  quick  setting  cement  and  the  surfaces  between 
the  end  blocks  and  the  center  block  were  greased  to  ensure  an  uniform  load 
distribution.  The  loading  was  displacement  controlled  by  using  two  Linear 
Variable  Differential  Transducer  (LVDT)  which  monitored  the  actual  deformations , 

In  order  to  measure  the  opening  and  sliding  of  the  crack  close  to  the 
initial  notch  tip,  two  clip  gages  specially  design  at  our  laboratory  were  used 
(Fig.  III. 4).  Three  specimens  of  each  group  (each  angle)  were  tested.  Clip 
gages  were  not  used  on  the  first  specimen  of  each  group  because  the  crack 
propagation  directions  were  not  known  and  because  the  investigators  were 
apprehensive  of  unstable  failure  damaging  the  gages.  The  clip  gages  were 
mounted  on  aluminum  mounts  as  shown  in  figure  III. 4.  The  concrete  surface  was 
first  polished  and  the  mounts  were  fixed  with  super  glue  adhesive.  Small  amount 
of  epoxy  was  used  on  the  side  of  the  mounts  to  provide  additional  strength.  The 
finite  dimension  of  the  aluminum  mounts  (about  1/4"  or  0.64  cm)  and  the  epoxy 
makes  it  to  average  the  crack  opening  and  sliding  over  a  finite  distance. 
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Hence,  the  clip  gages  effectively  measure  the  crack  opening  and  sliding  at  a 
distance  of  about  1/4"  from  the  notch  tip.  The  data  from  the  LVDTs,  the  Load 
Cell  and  the  clip  gages  were  automatically  monitored  by  a  four  channel  Nicolet 
4094  digital  oscilloscope. 

Ill . 5  Holographic  Interferometry 

The  experimental  focus  of  this  report  was  on  developing  an  adequate 
experimental  technique  for  studying  the  phenomena  described  above.  Since  the 
crack  propagation  in  the  experiments  conducted  were  at  different  places  on  the 
specimens  and  the  crack  path  was  inclined  and  not  known  a  priori  it  was 
necessary  to  have  a  technique  with  full  field  detection  capability.  It  is 
necessary  to  have  a  technique  with  the  highest  sensitivity  (fraction  of  a  micro¬ 
meter,  comparable  to  the  intercristal  space  in  mortar).  Due  to  the  nature  of 
the  material  and  the  size  of  the  specimens,  the  technique  should  be  noninvasive 
and  noncontacting  so  as  to  have  the  least  influence  on  the  specimen  tested. 
Holographic  Interferometry  (HI)  was  adopted  as  the  technique  that  best  satisfied 
all  the  criteria  described  above.  This  method  which  is  typically  used  for 
visual  inspection  of  cracking,  delamination  and  for  out  of  plane  displacement 
measurement  was  used  in  this  experiment  to  simultaneously  measuring  both  the  in 
plane  shear  and  tensile  crack  openings  by  using  an  innovative  technique. 

Two  specimen  of  each  group  were  tested  by  an  automated  holographic  system 
(Holographic  6000)  .  It  was  possible  to  make  a  large  number  of  holograms  by  this 
system  and  observe  the  propagation  of  crack  at  different  stages  of  loading1281 . 
Crack  propagation  starts  from  the  notch  tip  and  the  crack  curves  to  the 
direction  parallel  to  the  direction  of  applied  loading  (Fig.  III. 5).  The 
holograms  of  some  typical  specimens  are  shown  in  figure  III. 6.  Crack 
propagation  is  evident  as  discontinuities  in  the  fringe  pattern. 

A  specimen  of  each  group  was  tested  by  an  innovative  approach  based  on 
quantitative  measurement  techniques  using  holography1291 .  The  fringes  in  HI  are 
caused  by  the  displacements  on  the  object  surface.  Each  fringe  corresponds  to 
the  component  of  the  displacement  along  the  sensitivity  vector.  The  Sensitivity 
Vector  (S)  for  any  point  on  the  object  is  the  bisector  of  the  observation 
direction  (OQ)  and  the  illumination  direction  (IQ)  (Fig.  III. 7).  Hence,  a 
single  observation  gives  only  one  displacement  component.  In  order  to  find  the 
three  independent  displacement  components  of  a  specimen  it  is  necessary  to  use 
three  different  illumination  or  observation  directions.  Typically  this  is 
restricted  by  the  cost  of  additional  optics  and  the  trouble  involved  in  making 
three  times  as  many  holograms.  It  is  also  difficult  at  times  to  avoid  rigid 
body  motion  in  between  making  the  three  holograms  for  each  loading  stage.  In 
figure  III. 7  if  the  displacement  5  of  the  point  Q  has  a  component  along  S  equal 
to  QC,  the  path  length  of  the  beam  at  that  point  changes  by  the  amount  (QA-t-QB)  . 

QA  +QB  -  2  QC  cos 9  and  the  component  QC  -  5  cosip 

So  the  total  path  change  is  equal  to:  25  cos  S  cos  ip 
The  path  length  change  contributing  to  a  fringe  (dark  or  light  band)  is  A/2 
where  A  is  the  wave  length  of  the  He  Ne  laser  used  in  this  work  (.63 
micrometers).  Hence,  a  fringe  count  of  n  corresponds  to  a  displacement  of: 
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nA 

6 - 

4  cos 6  cos tJ> 

It  is  possible  to  count  fringe  to  an  accuracy  of  1/2  fringe  (from  the 
middle  of  a  dark  or  light  fringe  to  the  edge  of  the  same  fringe)  and  in  our 
experiment,  6  was  approximately  equal  to  25°.  The  accuracy  of  the  set-up  to 
measure  displacement  along  the  sensitivity  vector  is  limited  to: 

(1/2)  (0.63) 

accuracy  -  -  -  0.09  micrometers 

4  cos  30 

The  method  adopted  here  gets  around  these  draw  backs  without  any 
significant  sacrifice  in  accuracy.  This  is  achieved  by  using  a  single  4"  X  5" 
(10.2  X  12.7  cm)  holographic  plate  placed  very  close  to  the  specimen  surface 
(Fig.  III. 8).  This  makes  each  point  on  the  object  surface  subtend  a  large  angle 
on  the  holographic  plate.  Effectively,  it  is  now  possible  to  look  at  the  object 
through  four  corners  of  the  plate.  The  observed  fringe  pattern  is  different 
from  the  four  corners  due  to  the  differences  in  observation  directions.  By 
counting  the  fringe  patterns  from  the  four  corners  of  the  same  plate  we  have  a 
set  of  four  equation  involving  the  three  unknown  displacement  components  and  the 
four  different  sensitivity  vectors.  The  accuracy  of  this  set-up  is  demonstrated 
in  figure  III. 9.  If  the  plate  is  placed  close  to  the  specimen,  A  and  B  era  the 
two  points  on  the  plate  through  which  observation  of  the  fringes  are  made.  S' 
and  S"  are  the  two  sensitivity  vectors  corresponding  to  the  observation 
directions  A  and  B.  If  the  in  plane  displacement  of  point  Q  is  6 ,  its 
components  along  S'  and  S"  are  6  dot  S'  and  6  dot  S"  (dot  product).  6  is 
evaluated  from  the  difference  in  fringe  patterns  as  observed  from  A  and  B.  As 
in  figure  III. 8,  the  displacement  6  caused  a  path  length  change  of  26  cosd' 
cos and  26  cosd"  cos^"  respectively  to  the  observation  points  A  and  B 
respectively.  ij>'  and  \j>"  are  the  angles  between  the  in  plane  displacement  and 
the  sensitivity  vectors.  If  the  fringe  counts  at  points  A  and  E  were  nt  and  n2 
respectively,  then: 


4  6  cosd'  cost/)'  4  6  cosd"  cosx/>" 


the  difference  in  fringe  counts  from  A  and  B  is : 

(n1-n2)-4<S(cos0 '  cos^'-costf”  cos^")/A 

A  computer  program  was  developed  which  could  calculate  the  sensitivity 
vectors  when  the  coordinates  of  the  optics  were  input.  The  fringe  counts  were 
input  into  the  program  which  solved  the  four  equations  in  the  least  square  sense 
to  find  three  unknown  displacement  components.  The  in  plane  displacement 
components  (crack  opening  and  sliding  profiles)  were  then  evaluated  by  the 
program.  It  needs  to  be  mentioned  that  placing  the  holographic  plate  in  front 
of  the  specimen  makes  the  set-up  more  sensitive  to  the  out  of  plane  direction 
than  in  the  in  plane  direction.  It  is  necessary  to  place  the  plate  as  close  to 
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the  specimen  as  possible  to  increase  the  in  plane  sensitivity.  Proximity  of  the 
plate  to  the  specimen  is  restricted  by  other  considerations  such  as  the 
coherence  length  of  the  laser,  ratio  of  beam  intensities  (reference  beam:object 
beam),  etc.  Despite  the  typically  higher  out  of  plane  sensitivity  of  HI,  the 
out  of  plane  motions  computed  for  the  test  specimens  were  much  smaller  than  the 
crack  opening  and  sliding  displacements. 

111. 6  Experimental  Results 

The  loads  at  different  lengths  of  crack  propagation  were  noted  while  crack 
extensions  were  observed  by  HI.  A  FEM  study  was  done  to  find  to  find  the  SIFs 
and  crack  opening  and  sliding  displacements  corresponding  to  unit  load  at  1/2", 
1",  1.5"  and  2"  crack  extensions.  The  actual  experimentally  observed  crack  was 
analyzed  in  the  FEM  program  where  the  gradual  curvature  of  the  crack  was 
incorporated  by  a  series  of  linear  segments.  The  FEM  analysis  results  where 
multiplied  by  experimentally  observed  crack  propagation  loads  to  get  the  results 
plotted  in  figures  III. 10  and  III. 11.  The  critical  values  of  Kx  (or  KIe)  at 
which  cracks  propagated  is  seen  to  vary  widely.  There  is  also  a  large 
discrepancy  in  the  crack  profiles  shown  in  figure  III. 11.  While  for  the  36° 
specimen  the  crack  openings  were  mode  than  that  predicted  by  FEM  analysis,  those 
for  the  72°  specimen  were  less.  This  was  supposed  to  have  been  due  to  the 
traction  and  dilatancy  caused  by  intergranular  friction,  ligament  connection, 
etc.  It  was  presumed  that  while  the  aggregate  interlock  prevented  opening  of 
the  crack  for  the  72°  specimen,  sliding  induced  additional  opening  (dilatancy) 
in  the  36°  specimens. 

111. 7  Cohesive  Crack  Model 

In  order  to  account  for  the  phenomena  observed  in  figure  III. 10  and  III. 11 
an  analysis  based  on  the  model  described  in  part  I  was  used.  The  closing  and 
sliding  tractions  were  unknown  in  the  absence  of  any  adequate  model  and 
sufficient  experimental  evidence  of  mixed  mode  fracture  behavior  of  mortar 
cracks.  These  tractions  were  obtained  by  applying  unit  opening  and  closing 
forces  at  intermediate  points  on  the  crack  and  finding  the  displacement  profiles 
along  the  crack  due  to  that  unit  force.  If  four  points  along  the  crack  are 
considered,  this  results  in  a  8  by  8  matrix  M  where: 

[L]i*8[M]8*8  -  [D] 

D  is  obtained  from  the  difference  between  FEM  and  HI  displacements  (Fig. 
III.  11).  L  is  obtained  from  the  solution  of  the  equation  above.  So  L  was 
obtained  as  the  necessary  traction  forces  on  the  crack  faces  which  could  explain 
experimentally  observed  displacements.  The  values  of  KIc  were  corrected  for  the 
values  of  L  to  give  the  corrected  values  shown  in  figure  III.  12.  It  can  be 
observed  that  the  values  of  KIo  (corrected)  are  within  a  reasonably  narrow  band 
showing  the  possible  applicability  of  a  LEFM  based  criterion. 

1 1 1. 8  Conclusions 

1.  Holographic  interferometry  was  found  to  be  an  excellent  technique  to  study 
real  time  crack  propagation  with  high  sensitivity. 

2.  A  single  plate  HI  method  can  retain  the  accuracy  of  HI  while  measuring  in 
plane  displacements.  This  method  eliminates  lot  of  the  expenditure  and 
labor  involved  in  multiple  beam  HI  and  can  evaluate  full  field  3-D 
displacement  field  with  a  sensitivity  of  about  0.2  micrometers. 
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3. 


A  FEM  program  based  on  the  Maximum  Hoop  Stress  (MHS)  criteria  can  predict 
the  path  of  crack  propagation  in  mortar  under  compression- shear  loading 
with  reasonable  accuracy. 

4.  The  MHS  criterion  can  not  predict  crack  propagation  loads  due  to  the 
complexity  involving  interface  tractions.  The  effects  of  these  tractions 
have  to  be  accurately  modeled  and  considered  if  a  LEFM  based  crack 
propagation  criterion  is  to  be  postulated. 

5.  An  analysis  involving  the  interface  tractions  contributed  by  a  aggregate 
interlock,  ligament  connection  and  dilatation  holds  promise  as  a  means  of 
modeling  mixed  mode  fracture  behavior. 
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FIGURES 


Fig.  1.1.  Cohesive  Crack  Models.  (A)  Dugdale-Barrenblatt ,  (B)  Hillerborg, 

(C)  Jenq  and  Shah,  (D)  Cook,  et.  al. 

Fig.  1.2.  Center  Notched  Plate  Specimen. 

Fig.  1.3.  (A)  Loading  Arrangement,  (B)  Typical  Holographic  Interferogram. 

Fig.  1.4.  Comparison  of  LEFM  Crack  Profiles  with  Holographically  Measured 

Crack  Profiles. 

Fig.  1.5.  Comparison  of  Closing  Pressures. 

Fig.  1.6.  vs.  Crack  Length  for  Various  Load  Levels. 

Fig.  1.7.  Measured  Crack  Length  vs.  Predicted  Crack  Length  for  Kx-0. 

Fig.  1.8.  Crack  Profile  Comparison  for  K:  -  0. 

Fig.  1.9.  Closing  Pressures  Calculated  from  Known  Crack  Profiles  and  the 

Proposed  Bilinear  Closing  Pressure. 

Fig.  I. 10.  Crack  Profile  Comparison  fo*  Proposed  Closing  Pressure  and  K^O. 

Fig.  I. 11.  Kt  vs.  Crack  Length  for  Proposed  Closing  Pressure. 

Fig.  1.12.  Crack  Profile  Comparison  for  Proposed  Closing  Pressure  and  K^O. 
Fig.  II. 1.  Center  Notched  Mortar  Plate  Specimen. 

Fig.  II. 2.  Holographic  Set-up. 

Fig.  II.  3.  Image  Reconstruction  and  Acquisition  Set-up. 

Fig.  II. 4a.  Intensity  Distribution  Across  an  Interferogram. 

4b.  Isolation  of  Interferometric  Effect. 

4c .  Binary  Image . 

Fig.  II. 5.  Fringe  Count  by  Difference  in  Intensity  Value. 

Fig.  II. 6.  Holodiagram. 

Fig.  II. 7.  Holographically  Measured  Crack  Profiles. 

Fig.  II. 8.  Holographically  Measured  Strain  Fields. 

8a.  Crack  length  1.375  in;  load  -  1324  lb. 

8b.  Crack  length  1.875  in;  load  -  1461  lb. 

8c.  Crack  length  2.25  in;  load  -  1601  lb. 
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Fig.  II. 9.  Fracture  Process  Zone 

9a.  Crack  length  1.375  in;  load  -  1324  lb. 

9b.  Crack  length  1.875  in;  load  -  1461  lb. 

9c.  Crack  length  2.25  in;  load  -  1601  lb. 

Fig.  II. 10.  Difference  Between  Holographic  and  LEFM  Strain  Fields. 

10a.  Crack  length  1.375  in;  load  -  1324  lb. 

10b.  Crack  length  1.875  in;  load  -  1461  lb. 

10c.  Crack  length  2.25  in;  load  -  1601  lb. 

Fig.  II. 11.  Comparison  Between  Holographically  Obtained  and  Computed  Crack 
Profiles . 


Fig.  11.12.  Difference  Between  Holographic  and  Bilinear  Closing  Pressure 
Strain  Fields. 

12a.  Crack  length  1.375  in;  load  -  1324  lb. 

12b.  Crack  length  1.875  in;  load  -  1461  lb. 

12c.  Crack  length  2.25  in;  load  -  1601  lb. 

Fig.  III.l.  Infinite  Plate  and  Actual  Specimen. 

Fig.  III. 2.  Finite  Element  Mesh. 

Fig.  III. 3.  Comparison  of  Different  Solutions  for  Kx. 

Fig.  III. 4.  Clip  Gages  and  Mounts. 

Fig.  III. 5.  Tested  Specimens. 

Fig.  III. 6.  Holographic  Fringes. 

Fig.  III. 7.  Sensitivity  of  Holographic  Interferometry. 

Fig.  III. 8.  Four  Observation  from  One  Plate. 

Fig.  III. 9.  In  Plane  Displacement  Measurement. 

Fig.  III. 10.  Values  of  KIc  Before  Tractions. 

Fig.  III. 11.  FEM  vs.  Holographic  Crack  Profiles. 

Fig.  III. 12.  Corrected  Value  of  KIc. 
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Fig.  1.5.  Comparison  of  Closing  Pressures. 


Various  Load  Levels 
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Fig.  II. 2.  Holographic  Set-up. 
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Fig.  II. 3.  Image  Reconstruction  and  Acquisition  Set-up. 
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Fig.  II. 7.  Holographically  Measured  Crack  Profiles. 


Fig.  II. 8.  Holographically  Measured  Strain  Fields. 
8a.  Crack  length  1.375  in;  load  -  1324  lb. 

8b.  Crack  length  1.875  in;  load  -  1461  lb. 

8c.  Crack  length  2.25  in;  load  -  1601  lb. 
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Fig.  II. 9.  Fracture  Process  Zone 

9a.  Crack  length  1.375  in;  load  -  1324  lb. 
9b.  Crack  length  1.875  in;  load  -  1461  lb. 
9c.  Crack  length  2.25  in;  load  -  1601  lb. 


Fig.  II. 10.  Difference  Between  Holographic  and  LEFM  Strain  Fields. 

10a.  Crack  length  1  375  in;  load  -  1324  lb. 

10b.  Crack  length  1.875  in;  load  -  1461  lb. 

10c.  Crack  length  2.25  in;  load  -  1601  lb. 
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12c.  Crack  length  2.25  in;  load  -  1601  lb. 
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Fig.  III.l.  Infinite  Plate  and  Actual  Specimen. 
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Fig.  III. 5  Tested  Specimens 


Fig.  III. 6  Holographic  Fringes 
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Fig.  III. 11.  FEM  vs.  Holographic  Crack  Profiles. 
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